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C i Abstract Interaction between dark energy and dark matter is probed through 

deviation from the virial equilibrium for two relaxed clusters: A586 and A1689. 



The evaluation of the virial equilibrium is performed using realistic density pro- 
files. The virial ratios found for the more realistic density profiles are consistent 
with the absence of interaction. 
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i-i , 1 Introduction 

^ \ Searching for evidence and possible strategies for detection of dark energy and 

dark matter are among the most pressing issues of contemporary physics. Dark 
CSJ ' energy (DE) and dark matter (DM) are the fundamental building blocks of the 
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cosmological standard model based on general relativity. If on one hand, the 
interaction of the dark sector with the standard model such as neutrinos [T] 
and the Higgs field [2] lead to well defined experimental implications, such as 
for instance, a possible deviation from the Rutherford-Soddy radiative decay 
law ^ (see eg. Ref. [S] for general discussions), on another hand, interaction 
within the dark sector itself should not be excluded. Indeed, a recent study [B] 
suggests that a putative interaction might also be detectable using for instance 
gamma-ray bursts. The interaction between dark energy and dark matter can 
be introduced either by ad hoc arguments [7] or because dark energy and dark 
matter are unified in the context of some framework [51I51ITU]. Of course it is 
quite relevant to search for observational evidence of this interaction. 

In a previous work jllj we have analysed the effect of the interaction be- 
tween dark energy and dark matter on the virial equilibrium of clusters. Our 
method is based on the generalisation of the Layzer-Irvine equation, the cos- 
mic virial theorem equation, which was then applied to the Abell cluster A586. 
This cluster was chosen given that it is presumably in equilibrium and has not 
undergone interactions for quite a long time |12j . We have also argued that, 
based on the analysis of the bias parameter, that the dark sector interactions 
do signal a possible violation of the Equivalence Principle at cosmological scale 
[mil3) . This proposal has been extended for other clusters in Refs. [Hj, and 
as in our A586 analysis, evidence for DE-DM interaction was encountered. 

In this work we reexamine the cluster A586 and extend our analysis to the 
Abell cluster A1689 considering that mass distribution within the cluster is 
not constant. The additional cluster A1689 is chosen as it shares with A586 
the feature of showing a negligible amount of mergers, as inferred from its low 
X-ray substructure [13]. We consider four distinct density profiles: Navarro, 
Frenk and White (NEW) profile [H], isothermal sphere (ISO) profile fl71ll8[ 
milinilllllSa, Moore's (M) profile ^ and Einasto profile [24,25,26 . We also 
obtain results for up to four different techniques of evaluating the velocity 
dispersion (t„ of the haloes. We argue that the most relevant technique for 
determination of the velocity dispersion is through the X-ray temperature. 
Our analysis reveals that the use of more realistic mass profiles brings the 
virial ratio close to its canonical value in the absence of interaction. 

This work is organized as follows: in section [2] we describe the techniques 
used for fitting the data with the various density profiles; section [3] presents 
our analysis of the virial ratios for the previously used cluster A586; section |3] 
considers the A1689 cluster; in section [5] we interpret our results in terms of 
DE-DM interaction. Einally, in section |6] we present our conclusions. 
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2 Computations of the Virial ratio for a given density profile 

In this section we present our method and summarise the results of applying 
realistic density profiles to estimate the interaction between DE and DM in 
relaxed systems as discussed in Refs. [TT1IT5] . We will focus on two relaxed 
galaxy clusters: A586, which was already analysed in Ref. with a top-hap 
density profile, and A1689. 

We start by defining the basic quantities and deriving the necessary rela- 
tions between them. In particular we focus on obtaining various density profile 
fits for the haloes, and using the photometric coordinates of the galaxies, we 
compute the average distance between a galaxy and the centre of the cluster, 
(R) , at which point we are in condition to evaluate the various density profiles 
parameters and the corresponding errors for the cluster A586. For the cluster 
A1689, we use weak lensing surface density data. To compute the virial ratio 
Pk/pw we must also find the velocity dispersion ct„. From Refs. p^l28l[27l[29] 
we obtain the weak lensing and photometric data necessary to compute the 
different density profiles as well as values for (T„, as given by different methods; 
note that we approximate the velocity dispersion by its average value therefore 
considering it to be constant over the cluster ; we can then compute the ratio 
of kinetic energy to potential energy using different data sources (see Table [1] 
for A586 and Table [3] for A1689). 

In Ref [TT] we have used the simplest approach of a Top- Hat density profile 
and used the velocity dispersion estimated from weak lensing measurements. 
For the purpose of comparison, we present in Table [5] the outputs of such 
profile with the same velocity dispersion as for the other more realistic profiles. 
Although that analysis allows for claiming detection [11], a more accurate 
estimate of the density leads to an underestimation of the potential energy by 
placing more mass far from the centre than a realistic decreasing profile. This 
implies an overestimate of the magnitude of the virial ratio, as can be seen 
for central values in Table [21 This should be corrected with the use of more 
realistic density profiles. 

It is important to stress that the methods used vary slightly for the two 
clusters we are studying here. For the A586 we proceed in the same spirit as in 
Ref. [TT], estimating the density and shape parameters for the several density 
profiles from total 3D mass and galaxy coordinates. On the other hand, for the 
A1689 cluster we fit directly the profiles to the 2D surface density obtained 
from lensing measurements [28) . 



2.1 NFW profile 

Given its prominence in discussions about the density profile emerging from 
realistic N-body simulations in the context of the ylCDM model, the NFW 
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density profile 16 will be used. It reads 



p(r) 



^ 1 + - 



(1) 



where po and tq are the density and shape parameters respectively. Note that 
since p only depends on r, spherical symmetry is built into these computations 
from the very start. 

In the analysis of the cluster A1689, the 2D surface density from lensing 
measurements is used to perform a fit in order to obtain the estimates for po 
and tq. For the cluster A586 we proceed by computing the mass by integrating 
Eq. ID) over the volume: 
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which after considering Eq. ([2|) can be written as follows 
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We can then get the shape parameter tq for A586 by numerically inverting 
(R) . We then use it to define the NEW density parameter pQ from the observed 
mass M contained within the radius R by inverting Eq. ([2]): 
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The kinetic energy density is defined in terms of the average velocity dis- 
persion, (T^, to be 
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Assuming a constant average velocity dispersion, it reads 
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as used in Ref. [TT] 
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The potential energy is given by 
47r 
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(9) 



As discussed in Ref . [TT] , the existence of interaction between DE and DM 
is estimated by comparing the ratio pk/pw with the expected —1/2 value 
arising from the virial theorem. Taking the ratio of Eqs. ([7]) and © we find: 
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Throughout our analysis, the errors of the ratios are computed according to 
the following equation: 



A ( ^({oJ) ) = 
Pw 



\ Pie{o, 



f dpx/pw 



(11) 



where {oi\ are the parameters carrying measurement errors and are defined 
for each cluster as {oi} = {ay,M} for A586 and as {oi} = {a^,ro,po} for 
A1689. 



2.2 Isothermal profile 

Introduced first as the natural outcome of spherically symmetric DM self- 
gravitating infall P71I18I[T^I201I211[^ , the isothermal density profile is simpler 
than the NFW one: 

^° . (12) 



Here the fiducial radius rg and density poi or mass Mq = 47r/3/9o''oi are ar- 
bitrary as the profile is self-similar: there is no characteristic scale, so we can 
chose the mass and total radius of the halo as the fiducial values. 

Mo = M, (13) 
ro = R. (14) 

The mass is found by integrating Eq. (|T^ over the volume: 

/R 2 n 

-^-^dr = inrlRpo = Mq — . (15) 
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The mean radius (R) can then be defined with 

M {R) = An / -j^dr = 27rrlpoR\ 
° (^) 

which after considering Eq. (jlSp can be written as follows: 

The new kinetic energy density is now defined to be 
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Assuming a constant average velocity dispersion, it reads 
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again as in Ref. [llj . 

The potential energy is given by 
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Taking the ratio of Eqs. ^ and 1^ we find: 
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2.3 Moore's profile 

We consider now the Moore density profile that also arises as a suitable dark 
halo profile from N-body simulations [13]. It is believed to be quite accurate 
to describe galaxy size halo formation |30) : 
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where po and tq are the new mass density and shape parameters, respectively. 
The mass is found by integrating Eq. (|22p over the volume: 
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The mean radius (R) can then be defined with 



M {R) = inrlpo / ^ , , dr, (24) 
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and considering Eq. (j23p . it can be written as: 
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We can get the shape parameter ro numerically by inverting (R), which 
can be then used to define the Moore density parameter, po, with the observed 
mass, M, contained within the radius, i?, by inverting Eq. (123^ : 

3M 

Po 



8.r3ln(^l+(ii)^ 
Thus, the kinetic energy density is defined to be 
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Assuming a constant average velocity dispersion, it reads 
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still as in Ref. [TT]. 

The potential energy is given by 
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Taking the ratio of Eqs. 



and (HHI we find: 
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Table 1 Velocity dispersions from the various observations of A586 as given by Ref. 1121 . 



Method 


a (Km/s) 


X-ray Luminosity 


1015 ± 500 


X-ray Temperature 


1174 ± 130 


Weak lensing 


1243 ± 58 


Velocity distribution 


1161 ± 196 



2.4 Einasto's profile 



The Einasto density profile [Mll^l^ was originally used to describe the in- 
ternal density profiles of galaxies and has been proposed as a better fit model 
for ylCDM haloes EHEIES]: 



p = Pe exp 




Po exp 
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(31) 



where po — Pe^'^" is the central density and r_2 the radius at which the slope 
d\n p/d\nr = —2, the isothermal value. The radius r^. is defined such that it 
contains half the total mass and d„ is an integration boundary to ensure that. 
n gives the strength of the density fall. From Eq. (pij) the corresponding mass 
then reads: 



M = 47r 



Po exp 



-2n 
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r dr. 



(32) 



The mean radius {R) then becomes defined with 



M (i?) = 47r / Po exp 
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(33) 



We can get the shape parameter r_2 from observing the mean intergalactic 
distance and numerically solving Eq. ([55)) together with Eq. ([5^. 

We can then use it to compute the central density parameter po by making 
use of Eq. (l32|) . The kinetic and potential energy densities are computed from 
their definition 
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(34) 



through numerical integration. Note that it is also possible to perform the 
integrations analytically as done for the previous cases and find expressions 
relating the profile parameters and the virial ratio to the data M, (i?) and a. 
However, since these expressions for the Einasto profile are not particularly 
illuminating we will omit them. 
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Table 2 Virial ratio from the various observations of A586 obtained from data of Ref. |12| 
using different density profiles. 



Method 



X-ray Luminosity 
X-ray Temperature 
Weak lensing 
Velocity distribution 



Top Hat 



-0.516 ± 0.516 
-0.691 ± 0.190 
-0.774 ±0.145 
-0.676 ± 0.253 



NFW 



-0.408 ±0.407 
-0.545 ±0.150 
-0.611 ±0.115 
-0.533 ±0.200 





Moore 


Einasto n = 1 


Einasto n = 6 




-0.316 ±0.315 


-0.416 ±0.415 


-0.405 ± 0.404 




-0.423 ±0.116 


-0.556 ±0.153 


-0.542 ±0.149 




-0.474 ±0.089 


-0.624 ±0.117 


-0.608 ±0.114 




-0.413 ±0.155 


-0.544 ± 0.204 


-0.530 ±0.199 



Isothermal 



-0.353 ±0.352 
-0.472 ±0.130 
-0.529 ±0.099 
-0.461 ±0.173 



3 Analysis of A586 cluster 

3.1 Data Analysis 

In the analysis of A586 we used the same data used in the original analysis 
performed in FTI], that is: 

— member galaxy coordinates 

— total mass inside a radius R, obtained from weak lensing measurementes 

M = (4.3 ± 0.7) X 1O"M0 R = 422 Kpc (35) 

— velocity dispersion from different sources, table [T] 

Using the coordinates of the galaxies that compose A586 we can compute the 
average distance between a galaxy and the centre of the cluster, (i?). We start 
by determining the centre of the cluster. This is done by computing the average 
declination and right ascension from the coordinates of the 31 galaxies in the 
cluster. 

The distance of a galaxy i, with coordinates {ai,Si), to the centre of the 
cluster (ac,<5c) is given by: 

= 2cP[l — cos{ai — ac)cos{Sc)cos{Si) — sin{5c)sin{5i)]. (36) 

Using Eq. (1361) to compute for all galaxies in the reduced sample and then 
taking the average: 

{R)=22i.QKpc. (37) 

We can then proceed to compute the shape parameters of the various 
density profiles as described in the previous section and through them compute 
the virial ratio and the interaction between DE and DM. 
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Table 3 Summary of A1689 available data 



Parameter 



Value 



Reference 



po (10 ''^h''gr/crw 
ro{h~^kpc) 
kT^KeV) 

(^dyn{Km/s) 

Ray„(h-'^Mpc) 



7) 9.6 ± 1.8 

175 ± 18 




2.26 



M 

[28] 

[23 

[29] 

'29' 



Table 4 NFW fit to A1689 lensing data. 



Parameter Value Error 



po (10-^" /i^gr/cm^) 4.513 0.330 
ro{h-^kpc) 271.6 12.5 

■R'^ 0.999 



3.2 Comments on the results for A586 

We present our results in table [2] One sees that the analysis of this cluster 
with a top hat profile is in agreement with the findings of Ref. [TT]. We ob- 
serve consistently that using weak lensing velocity dispersion (as performed in 
Ref. [TT]) yields a higher virial ratio. We stress that the weak lensing velocity 
dispersion is not the most adequate data source as it introduces correlations 
between mass and velocity estimation. It is crucial to avoid these correlations 
as the aim of this work is to detect deviations to the virial equilibrium and 
interpret these as an effect due to DE-DM interaction. 

4 Analysis of the A1689 cluster 

In table [3] we gather the data available for A1689. As already mentioned, for 
this cluster we find the density and shape parameters for the different profiles 
by directly fitting them to the 2D surface density obtained from weak and 
strong lensing. The relation between the 2D surface density k{R) and the 
total 3 dimensional mass density p(r) is given by 



where Scrit is the critical density for lensing, which depends on the back- 
ground cosmology and also on the lens and source sample. For the data we 
are considering, Scrit — 1.0122 /ig/cm^ [28], for Hq = lOOhkm/s/Mpc and 
observationally h ~ 0.7. Then, by letting p(r) be NFW, isothermal sphere, 
Moore, or Einasto, we estimate the parameters of these four density profiles. 

In the ensuing analysis, we will consider two methods for estimating veloc- 
ity dispersion: galaxy dynamics and X-ray temperature, as shown in table |31 
While the velocity dispersion from the galaxy dynamics is a direct measure- 
ment, converting the X-ray temperature into a velocity dispersion involves 
assumptions that must be explained. Defining j3 ~ a'^{kT / ^rrip)^^ , where ^ 




(38) 
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Table 5 Moore fit to A1689 lensing data. 



Parameter 


Value 


Error 


po (10-^''h^gr/cm^) 


1.064 


0.525 


ro{h~^kpc) 


357 


102 


7^2 


0.973 





Table 6 Einasto fit to A1689 lensing data. 





n = 


1 


n 


= 6 


Parameter 


Value 


Error 


Value 


Error 


po (10-''^ h'-'gr/cm-^) 


46.49 


5.165 


66005 


7358 


ro{h~^kpc) 


122.5 


8.47 


467.3 


31.99 


■R? 


0.9834 




0.998 





and nip are the reduced nuclear mass and the proton mass, allows for com- 
puting knowing kT once (3 is specified. To first approximation one can 
assume the hypothesis of density energy equipartition, which corresponds to 
setting (3 = 1. For this choice we find = (1232 ± 27) km/ s. Other values 
for (3 might be used, in particular see e.g. Ref. Another issue related to 
the velocity dispersion is the region in which it is measured. In the case of 
(jy coming from galaxy dynamics, the measurements are taken in a spherical 
region whose radius is the distance to the centre of the outermost galaxy. For 
the current cluster, this is found to be 2.26h~'^Mpc. In the case of X-ray tem- 
perature measurements, these are taken up to a radius ^2500, i-e. the radius at 
which the local density is 2500 times the critical density around the redshift 
of the cluster. Assuming the ACDM cosmology, we find this to be of order 
300/i~^A:pc. Notice that the exact value of r25oo depends on the density profile 
as it is assumed, nonetheless there is an order of magnitude difference relative 

to Rdyn • 



4.1 NFW profile 

In Ref. [28| a joint X-ray and lensing analysis is performed and the data has 
been fitted to an NFW profile, the results of which are displayed in table |31 
Nonetheless we perform our own fit to the lensing data, the results of which 
are shown in table |4l 

Since the method used for computing px / Pw for A1689 is slightly different 
from the one used for A586, the value and error of the ratio will now be given 
in terms of the fitted parameters po and tq instead of M and (i?) . The kinetic 
to potential energy density ratio is now given by 





[(' + 


^1 

ra J 


In 
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(-S) 




SnGpor"^ 
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- In 
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2 



(39) 



and the error is found by applying Ec^. (IllD to the previous equation. 
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Table 7 Virial ratio from the various observations of A1689 given by Refs. I28II27II29] . 



Method 


NFW 


NFW M. 


Isothermal 




X-ray Temperature 
Velocity distribution 


-0.518 ± 0.060 
-0.408 ± 0.099 


-0.538 ±0.144 
-0.485 ±0.173 


-0.715 ±0.067 
-0.647 ±0.146 







Moore 


Einasto n = 1 


Einasto n = 6 




-0.549 ±0.391 


-0.527 ±0.093 


-0.520 ±0.086 




-0.479 ±0.389 


-0.473 ±0.130 


-0.353 ±0.098 



For completeness we compute the ratio px/pw using our best fit and the 
fit found in Ref. [28] . The results are displayed in table [71 



4.2 Isothermal sphere density profile 

Following the procedure described above, we fit the projected mass data ob- 
tained from lensing measurements to the isothermal mass profile given by Eq. 
(IT2]) . One should notice that the isothermal profile, unlike the NFW, Moore 



and Einasto profiles, has only one free parameter, namely por^. The fit to the 
data yields por^ = (3.987 ± 0.333) x IQ-^i with TZ^ = 0.845. As indicated by 
the low value of TZ^ the isothermal profile is not very suitable to describe the 
mass distribution in the cluster A1689. This can also be seen in Figure[T] The 
energy density ratio is given in terms of the fitted parameter and the velocity 
dispersion by: 

---l?^^ (40) 
pw 87r Gr^po 

with an error given by Eq. ([TT|) . Note that Eq. pOI) is independent of the 
region where the velocity dispersion is measured, R. This is a characteristic 
feature of this profile. In Table [7| we display the results found by applying the 
isothermal profile to the study of A1689. 



4.3 Moore's profile 

Repeating the above procedure for the Moore's profile with data depicted in 
Table [S] we get the results shown in Table [7| from which we can see that no 
signal of DE-DM interaction can be seen as no deviation from the virial ratio 
value —1/2 is unambiguously seen. 



4.4 Einasto's profile 



Repeating the above procedure for the Einasto profile with data depicted in 
Table [51 we obtain the results shown in Table [7| 
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2 3 4 5 6 7 



Log[r (A ' Kpc)] 

Fig. 1 Fits to 2D surface density for A1689. 

5 DE-DM Interaction 

In this section we relate the virial ratios computed above to the interaction 
parameter between DM and DE foUowing Ref. [11]. We start by briefly re- 
viewing how the interaction is parametrised and how this parameter is related 
to the virial ratio in relaxed structures. In Ref. [TT] the interaction was anal- 
ysed within the context of two distinct models of interacting DE: coupled 
quintessence [7] and the generalized Chaplygin Gas (GCG) [5]. It was shown 
that for (jJde = — 1, the coupling parameter of the interacting quintessence, 
^ (cf. below), is related to the a parameter of the GCG equation of state 
p = —Ajp"-, by the scaling parameter (defined below) r\ — 3(1 + a). Keep- 
ing this mapping between the two distinct models in mind, we review how 
to relate virial ratio the interaction parameter in the context of the coupled 
quintessence model. 

The conservation equations for DM and DE read 

PDA/ + iHpDM = ^HpDAU (41) 



Pde + 3HpDE{i + i^DE) = -^HpDM- (42) 

It is assumed that there is a scaling behavior between the DE and DM 
energy densities, 

PD^ = (43) 

RDM i^DM 

where rj describes the scaling and is related to the coupling ^ by [TT] 

£ = (44) 

l + {{2DMo/nDEo)a-'^ 

Following the derivation of the Layzer-Irvine equation in Ref. |34| , noting 
that the scaling of the matter kinetic energy with the scale factor is unchanged 
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by the presence of interaction with DE, pK oc a and that the interaction 
will modify the scaling of the potential energy with a to pw oc a^~^ one finds 

m 

PDM + H{2pK + Pw)=^PwH. (45) 

This is the modified Layzer-Irvine equation. Assuming that locally p 0, we 
find that the virial theorem is modified to 

---hi (46) 

Pw ^ ^ 

where ^ is the interaction strength as defined in Eqs. (|4T|) and (j42p . 

We are therefore mapping deviations from virial equilibrium to the exis- 
tence of DE-DM interaction. If the data on a particular cluster yields px / Pw > 
— 1/2=>^>0 and the energy transfer occurs from DE to DM. Conversely 
PkI Pw <— 1/2=>^<0 and the energy transfer occurs from DM to DE. 

As already stressed in the previous section, the results in tables [2] and [7] do 
not point to a preferred direction of the deviation of the ratio from its canonical 
value of 1/2. Therefore there is no preferred sign for ^ and the results are largely 
consistent with ^ = 0. This contrasts with the results found in pUfTB] where 
the analysis of the A586 cluster with a top hat profile and weak lensing inferred 
velocity dispersion pointed towards ^ < 0. In fact, this suggests that changing 
the density profile and considering the related methodological difhculties do 
not allow for a definite conclusion on the virial ratio. We believe that the 
current analysis represents a step forward relative to the original proposal as 
it deals with more realistic density profiles and uses uncorrelated data sources 
in the estimation of the velocity dispersion. However, in order to distinguish 
between the different physical scenarios ^ > 0, ^ < and ^ = it is crucial to 
have a better understanding both of the mass distribution within the clusters 
and a more importantly a more accurate knowledge of their dynamical state. 
This will allow to decrease the error bars and detect putative deviations of the 
ratio from —0.5. 

6 Conclusions 

In this work we have analyzed the virial equilibrium of clusters A586 and A1689 
using four density profiles, namely NEW, isothermal spheres and Moore's 
and Einasto's profiles. The method employed represents an evolution from 
the original proposal as it deals with more realistic mass distributions. For 
the A586 cluster we have only found evidence of deviation of the virial ra- 
tio PkI Pw = for the over-simplistic Top-Hat density profile, in accor- 
dance with [11]. In what concerns cluster A1689, we have not encountered a 
convincing evidence of this interaction except when using the uni-parametric 
isothermal sphere distribution. This could reflect the claims on A1689 that it 
is not a relaxed cluster, in particular, that it exhibits signs of triaxiality '36j. 
Despite the fact that search for interaction with more realistic mass profiles 
has not returned a clear signal in either cluster, one should stress that the 
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encountered error bars are still too large at present. In order to study the 
interaction between DE and DM via deviations from the virial equilibrium it 
is crucial to have a better understanding of the cluster mass distribution and 
of their dynamical state. We believe N-body simulations have a crucial role to 
play in this problem. 
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